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^Metodo per filtrare una sequenza di immagini 
digital!" 
DESCRIZIONE 

La presente invenzione riguarda 1 ' elaborazione di 
iinmagini digitali ed in particolare si riferisce ad un 
metodo per filtrare il rumore in una sequenza di inunagini 
digitali . 

Le inunagini digitali sono attualmente utilizzate in 
molteplici applicazioni quali, ad esempio, quelle relative 
ai dispositivi d' acquisizione tradizionali, come le 
foto/video caitiere digitali. E ' previsto inoltre un sempre 
piu diffuso utilizzo delle inunagini digitali in 
dispositivi di nuova generazione, ad esempio, nei 
terminali mobili di comunicazione multimediale. 

Esistono numerosi dispositivi o applicazioni che 
utilizzano iinmagini digitali in sequenza, cioe acquisite a 
breve intervallo temporale l'una dall'altra e 
rappresentanti approssimativarnente una stessa scena reale. 

La velocita di acquisizione della sequenza, ossia il 
numero di immagini acquisite in un dato intervallo di 
tempo, e variabile in base alle specif iche applicazioni; 
ad esempio, e molto elevata nelle videocamere digitali 
(circa 25 immagini al secondo) ed e invece piu bassa in 
terminali mobili di comunicazione (circa 15 immagini al 
secondo) che acquisiscono le immagini digitali e le 



trasmettono in tempo reale verso un terminale remote 

E' noto che i dispositivi d' acquisizione di iinmagini 
digitali, specialmente se includenti sensori di tipo CMOS , 
introducono intrinsecamente del rumore nelle immagini 
acquisite. 

Nelle sequenze immagini digitali il rumore, oltre a 
degradare la qualita delle immagini, riduce anche 
l'efficienza di codif ica/compressione . Comunemente, 
infatti, le sequenze di immagini acquisite sono destinate 
ad essere codif icate/compresse con tecniche di 
codif ica/compressione che operano, ad esempiq, in accordo 
alio standard MPEG o alio standard H263 e che sono ormai 
diffusamente utilizzate nella maggior parte dei 
dispositivi attualmente in commercio. 

In presenza di rumore, l'efficienza di codif ica / 
compressione e ridotta perche il rumore introdotto e, 
tipicamente, in forma di fluttuazioni casuali che riducono 
la ridondanza sia all'interno di un'immagine sia fra 
immagini temporalmente vicine. 

Esistono numerose tecniche di filtraggio, per ridurre 
o eliminare il rumore presente in una sequenza di 
immagini . 

Sono stati numerosi i tentativi di sviluppare 
tecniche efficienti di riduzione del rumore di una 
sequenza utilizzando differenti specif ici tipi di filtri. 



Ad esempio, sono noti filtri digitali di tipo passa-basso, 
filtri mediani, filtri spaziali adattativi e filtri 
temporali ricorsivi con o senza compensazione del 
movimento . 

Sono, inoltre, note tecniche che, per migliorare 
l'efficienza di riduzione del rumore in sequenze di 
immagini, utilizzano metodi ibridi che combinano un 
filtraggio digitale di tipo spaziale ad un filtraggio di 
tipo temporale. 

Le tecniche note di riduzione del rumore in sequenze 
di immagini, seppur soddisfacenti sotto diversi aspetti, 
presentano numerosi inconvenienti e problemi, legati ad 
esempio a prestazioni non adeguate, a onerosita 
d'elaborazione o ad un' implement azi one eccessivamente 
costosa a tal punto da poter essere dif f icilmente 
impiegate in dispositivi d' acquisizione portatili di tipo 
commerciale. 

Scopo della presente invenzione e di mettere a 
disposizione un metodo di riduzione del rumore, in una 
sequenza di immagini, che non presenti gli inconvenienti 
dei metodi della tecnica nota. Questo scopo si ottiene con 
un metodo per filtrare una sequenza di immagini digitali 
in formato CFA cosl come descritto nelle rivendicazioni da 
1 a 12. 

Costituiscono oggetto della presente invenzione anche 
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un filtro cosl come descritto nella rivendicazione 13 ed 
un dispositivo di acquisizione cosl come descritto nella 
rivendicazione 14 . 

Ulteriori caratteristiche e i vantaggi 

5 dell' invenzione saranno meglio compresi dalla seguente 
descrizione dettagliata di una sua forma d'esecuzione 
preferita, data a titolo indicativo e pertanto in nessun 
modo limitativo in relazione alle figure allegate, in cui : 
la figura 1 mostra lo schema a blocchi 
10 esemplif icativo di un possibile dispositivo d' acquisizione 
che implementa un metodo d' elaborazione secondo la 
presente invenzione, 

la figura 2 mostra lo schema di disposizione 
degli elementi filtranti di un sensore di tipo Bayer 
15 impiegabile nel dispositivo di figura 1, 

la figura 3 rappresenta schematicamente la 
successione di fasi di un particolare metodo 
d' elaborazione secondo la presente invenzione, 

la figura 4 mostra una maschera di selezione 
20 per selezionare pixel di colore verde impiegabile nel 
metodo d' elaborazione secondo la presente invenzione, 

la figura 5 mostra due maschere di selezione per 
selezionare pixel di colore rosso e pixel di colore blu 
impiegabili nel metodo d' elaborazione secondo la presente 
25 invenzione, 
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la figura 6 mostra in maggior dettaglio una 
delle fasi della successione rappresentata in figura 3, 

la figura 7 mostra un esempio di scelta di pixel 
in accordo ad una selezione di tipo DRT , 
5 la figura 8 e un grafico che illustra risultati 

sperimentali , e 

la figura 9, riporta schematicamente due 
arch.it etture di filtraggio, la prima in accordo alia 
presente invenzione, la seconda di tipo convenzionale . 
10 L' esempio di realizzazione preferita della presente 

invenzione si riferisce ad un dispositivo portatile atto 
ad acquisire sequenze di immagini digitali per 
applicazioni video, ed in particolare si riferisce al 
filtraggio del rumore di una sequenza di immagini 
15 acquisita con una video camera digitale. 

A tal proposito, si osservi che gli insegnamenti 
della presente invenzione possono essere estesi anche ad 
applicazioni differenti da quella cui si fa esplicitamente 
riferimento nella seguente descrizione, ad esempio 
20 all ' acquisizione di sequenze di immagini in terminali 
mobili multimediali di comunicazione di nuova generazione. 

In figura 1 e mostrata, in maniera molto schematica, 
mediante blocchi funzionali una videocamera digitale 1. La 
videocamera 1 include un blocco 2 di acquisizione, 
25 comprendente un sensore ottico 3 . 



II sensore 3, ad esernpio di tipo CCD (Charge Coupled 
Device) o CMOS (Complemetary Metal Oxide Semiconductor) , e 
un circuito integrato comprendente una matrice di celle 
fotosensibili , ciascuna delle quali serve a generare un 
segnale elettrico proporzionale alia quantita di luce che 
la investe durante l'intervallo di acquisizione . Ciascuna 
cella f otosensibile del sensore, a cui comunemente ci 
riferisce con il termine di pixel, corrisponde ad un 
rispettivo pixel nell' immagine digitale. 

In una realizzazione preferita il sensore 2 comprende 
un f iltro ottico CFA (Colour Filter Array) , per esernpio a 
matrice di Bayer. 

Com'e noto ad un esperto del settore, in un sensore 
con f iltro CFA per acquisire un pixel dell' irnmagine e 
disponibile una sola cella f otosensibile . Il sensore e 
ricoperto da un filtro ottico costituito da una matrice 
(ad esernpio, di Bayer) d'elementi filtranti, ciascuno dei 
quali e associato ad una cella f otosensibile . Ogni 
elemento filtrante trasmette (assorbendone una porzione 
minima) alia cella f otosensibile ad esso associata la 
radiazione luminosa corrispondente alia lunghezza d'onda 
della sola luce rossa, della sola luce verde o della sola 
luce blu, in modo da rivelare per ciascun pixel soltanto 
una componente. Lo schema di disposizione degli element i 
filtranti in un filtro di Bayer e rappresentato in figura 



2, dove con R,G,B, sono indicati elementi filtranti 

rispettivamente di colore rosso, verde, blu. 

La videocamera 1 include inoltre un blocco di 

conversione analogico/digitale A/D, indicato con 4, per 
5 tradurre il segnale elettrico generato da ogni cella 

f otosensibile in un valore digitale su un prefissato 

numero di bit (generalmente 8, 10 o 12 bit) . Si puo 

assumere, solo a titolo d'esempio e senza per questo 

introdurre alcuna limitazione, che nella presente 
10 invenzione il convertitore A/D 4 sia tale da codificare i 

segnali analogici in ingresso con valori digitali 

rappresentati con otto bit. 

In uscita dal blocco A/D 4, I 7 immagine digitale 6 in 

formato CFA (Colour Filter Array) , cioe e costituita da 
15 una sola componente cromatica (R,G o B) per pixel. Per 

questo motive, un solo valore digitale su un byte e 

associato ad ogni pixel . 

Un blocco di filtraggio 5, e tale da filtrare il 

rumore operando direttamente sulle immagini digitali della 
20 sequenza in formato CFA, in questo esempio, in formato di 

Bayer, producendo in uscita, per ogni immagine CFA 

rumorosa in ingresso, una corrispondente immagine CFA con 

rumore ridotto. 

Un blocco di pre-processamento PrePro 6, attivo prima 
25 e durante tutta la fase di acquisizione, e tale da 
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interagire con il blocco di acquisizione 2 per estrarre 
dalle immagini CFA filtrate alcuni parametri utili ad 
eseguire funzioni automatiche di controllo: funzioni di 
auto-focus, auto-esposizione, correzione dei difetti del 
5 sensore e bilanciamento del bianco. 

Un blocco 7 denominato IGP (Image Generation 
Pipeline) e preposto alio svolgimento di una fase di 
elaborazione per produrre, a part ire da una immagine 
digitale CFA, un' immagine digitale completa, ad esempio in 

10 formato YCrCb, in cui ad ogni pixel sono associati tre 
valori digitali (in totale, 24 bit) , corrispondenti ad una 
component e di luminanza Y, e a due component i di 
crominanza Cr e Cb. Questa trasf ormazione, nota con il 
termine di interpolazione del colore, comporta un 

15 passaggio da una rappresentazione dell ' immagine su un solo 
piano (piano di Bayer), contenente 1 ' inf ormazione relativa 
a diverse component i cromatiche, ad una rappresentazione 
su tre piani. 

Comunemente nelle fotocamere digitali, il blocco. IGP 
20 e realizzato attraverso un processore dedicato (processore 
CFA) , ad esempio in tecnologia VLSI (Very Large Scale 
Integration) . 

Preferibilmente, in questo esempio il blocco 7 IGP, 
oltre all' interpolazione, e anche tale da eseguire diverse 
25 funzioni, quali ad esempio 1 7 applicazione di effetti 



9 

speciali, la correzione gamma, lo scalamento, la 
stabilizzazione ed altre funzioni in genere variabili in 
numero e tipo da costruttore a costruttore. 

Seguono un blocco 8 di compressione/codif ica, in 
5 questo esempio di tipo MPEG (ma puo essere di altri tipi 
come, ad esempio, H263), ed una unita di memoria 9. 

Nella videocamera 1, durante la ripresa di una 
sequenza video, le immagini della sequenza sono acquisite 
consecutivamente attraverso il blocco di acquisizione 2, 
10 preferibilmente, a breve distanza temporale l'una 
dall'altra. Ad esempio, lo standard MPEG-4 richiede che 
siano acquisite quindici immagini al secondo. 

In seguito indicheremo con Img lf Img 2 , Img 3/ Img n - le 
Img n , Img n + lf „. le immagini acquisite in sequenza: Irngx 
15 rappresenta la prima immagine della sequenza acquisita in 
ordine di tempo, Tmg 2 rappresenta la seconda immagine e 
cosi via. 

Ciascuna immagine, dopo 1 ' acquisizione, viene 
elaborata dai blocchi successivi, in modo che anche nelle 
20 successive fasi di elaborazione le immagini siano 
elaborate in sequenza nell' ordine temporale di 
acquisizione. 

Le immagini della sequenza, dopo 1 ' acquisizione, sono 
convertite in valori digitali dal convertitore A/D 4. 
25 Le immagini digitali in formato CFA sono poi inviate 
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in ingresso al blocco di filtraggio 5 CFA NF per essere 
elaborate in accordo al metodo di filtraggio del rumore 
secondo la presente invenzione. Il blocco di filtraggio 
produce in uscita una sequenza di immagini CFA filtrate, 
5 rispettivamente f_lmg lf f_Img 2 , f_Img 3 , f_lmg n - 1 , f_lmg n , 
f_Img n+1 , aventi rumore ridotto. 

Le immagini CFA filtrate sono poi elaborate dal 
blocco di pre-processamento 6. 

Ciascuna immagine in formato CFA in uscita del blocco 
10 di pre-processamento G e inviata al blocco di elaborazione 
IGP 7 . In tale blocco le immagini sono sottoposte ad una 
fase di interpolazione del colore e trasformate dunque in 
immagini complete, ad esempio in formato YCrCb. 

La fase d' interpolazione del colore pud essere, ad 
15 esempio, effettuata con metodologie note al tecnico del 
ramo e, pertanto, risulta evidente dalla precedente 
descrizione. 

Successivamente le immagini, sono inviate al blocco 8 
*MPEG Encoder' 7 che produce in uscita una sequenza, o 
20 flusso, di immagini codif icate/compresse secondo una 
codifica di tipo MPEG. 

II flusso di immagini compresse "MPEG Stream" pud 
essere registrato in un'unita di memoria 9, oppure inviato 
ad una perif erica esterna, non mostrata in figura. 
25 In una realizzazione preferita, il metodo 
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d' elaborazione filtra la sequenza di immagini CFA, in 
quest 7 esempio di Bayer, immagine per immagine. 

Le immagini CFA sono filtrate pixel per pixel, con un 
ordine di scansione sistematica dei pixel da sinistra 
5 verso destra e dall'alto verso il basso. In particolare, 
per ogni pixel p n (x,y) di un' immagine Img n e calcolato un 
rispettivo pixel omologo filtrato f jp n (x,y) di una 
corrispondente immagine filtrata f_Img n . 

Per elaborare il pixel p n (x,y) dell ' immagine lmg n e 
10 calcolare il pixel filtrato f _p R (x,y) , il metodo della 
presente invenzione vantaggiosamente utilizza anche pixel 
filtrati appartenenti ad una immagine della sequenza 
precedentemente filtrata f_Img n -z . Piu in particolare 
utilizza 1' immagine f_Img n . 2 risultante dal filtraggio 
15 dell ' immagine Imgn-x che nella sequenza precede 
temporalmente 1 ' immagine da f iltrare Img n . 

Tipicamente, dunque, per eseguire il processo di 
filtraggio sono sufficient! tre buf fer d' immagine : due 
buffer d'ingresso contengono rispettivamente 1' immagine 
20 Img n da filtrare (immagine corrente) e 1' immagine 
precedentemente filtrata f_Img n -i, e un buffer d'uscita 
contiene 1' immagine corrente filtrata f_Img n . 

In figura 3 e rappresentata schematicamente la 
successione di fasi di un metodo di elaborazione 2 0 per la 
25 riduzione del rumore in accordo alia presente invenzione. 
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In particolare, in figura 3 sono mostrate le fasi che 
consentono di ottenere a partire dal pixel p n (x,y) 
dell' immagine Img n il rispettivo pixel omologo filtrato 
f _J>n(XfY) della corrispondente immagine filtrata £_Img n . 

Una prima fase di selezione SW_SEL 21, dato il pixel 
in ingresso p n (x,y) da filtrare, seleziona un primo 
insieme di pixel SW n (x,y) comprendente il pixel p n (x f y) 
stesso ed una pluralita di pixel, appartenenti 
all ' immagine Img n , ad esso limitrofi. Piu in dettaglio, 
sono selezionati pixel limitrofi aventi associato lo 
stesso colore (R,G o B) del pixel da filtrare. 

In una realizzazione preferita, la selezione e 
effettuata utilizzando maschere, o matrici, di selezione 
SM_G, SMR, SMB differenti in base al colore del pixel da 
filtrare p n (x,y) , tutte di dimensione 5x5 , ad esempio come 
quelle rappresentate nelle figure 4 e 5. 

In figura 4 e mostrata una maschera di selezione SM__G 
per i pixel verdi (G) , in accordo ad una forma 
realizzativa dell' invenzione. Nella fase di selezione 21 
la maschera SM__G e allineata all' immagine da filtrare in 
modo che G 0 sia in corrispondenza del pixel verde p n (x,y) 
da filtrare. In tal modo, la maschera SM_G seleziona un 
primo insieme SW n (x,y) di pixel comprendente il pixel 
verde p n (x,y) in corrispondenza di G 0/ ed otto pixel verdi 
limitrofi disposti rispettivamente in posizioni 
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corrispondenti ai pixel Gi,...,G 8 della maschera SM__G 
rappresentata in figura. Tale insieme definisce una 
finestra operativa spaziale SW n (x,y) per il pixel verde da 
f iltrare. 

Analogamente, con riferimento alia figura 5, sono 
most rate le maschere di selezione SM_R ed SM_B, 
rispettivamente utilizzate se il pixel da filtrare e di 
colore rosso o di colore blu. 

La maschera SM_R, per i pixel rossi , e tale da 
selezionare un primo insieme SW n (x,y) di pixel 
comprendente il pixel rosso da filtrare p n (x,y) in 
corrispondenza di R 0 , ed otto pixel rossi limitrofi 
disposti rispettivamente in posizioni corrispondenti ai 
pixel Ri,...,R 8 della maschera. 

Si noti che, in questa particolare forma 
realizzativa, la maschera di selezione SM_B per i pixel 
blu e identica alia maschera di selezione per SM_R per i 
pixel rossi. Questa scelta e possibile grazie al 
particolare schema di disposizione degli element i 
filtranti in un sensore di tipo Bayer. 

Questo vantaggiosamente consente, durante la fase di 
selezione 21 SW_jSEL, di dover discriminare solo fra due 
possibili casi, cioe distinguere se il pixel p n (x,Y) da 
filtrare e verde o meno. 

Ritornando al diagramma di figura 3, una seconda fase 
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di selezione 22 TW__SEL, seleziona un secondo insieme di 
pixel TW n (x, y) , comprendente pixel appartenenti 
all ' immagine precedentemente filtrata f_Img n -i disposti in 
posizioni corrispondenti, cioe omologhi, ai pixel del 
5 primo insieme TW n (x,y) . 

Vantaggiosamente, a tal scopo, possono essere 
utilizzate le matrici di selezione sopra descritte 
applicate, questa volta, all' immagine precedentemente 
filtrata f_Jjng n _i. 
10 L' insieme di pixel cosi ottenuto definisce una 

finestra operativa temporale TW n (x,y) per il pixel da 
filtrare. 

Le finestre operative temporali e spaziali 
rappresentano 1' insieme dei pixel che interverranno nelle 
15 successive fasi del processo di filtraggio del pixel 
p n (x,y) . 

Come d not o ad un esperto.nel settore, nel filtraggio 
digitale di sequenze di immagini, un pixel filtrato pud 
essere ottenuto combinando in modo opportuno un certo 
20 numero di pixel ad esso adiacenti spazialmente (filtraggio 
spaziale) , temporalmente (filtraggio temporale) o spazio/ 
temporalmente (filtraggio spazio/temporale) . 

In particolare, come sara descritto meglio in 
seguito, il metodo della presente invenzione decide, pixel 
25 per pixel, se utilizzare un filtraggio di tipo 
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esclusivamente spaziale oppure utilizzare un filtraggio di 
tipo spazio/temporale . Vantaggiosamente la decisione sul 
tipo di filtraggio e legata all'entita del movimento fra 
immagini successive della sequenza consentendo di evitare 
5 la compensazione del movimento, di per se molto onerosa 
computazionalmente . 

Nel caso in cui sia utilizzato il filtraggio 
spaziale, in questo interverranno i pixel della finestra 
operativa spaziale SW n (x,y) , altrimenti saranno utilizzati 
10 pixel appartenenti ad entrambe le finestre. 

Una volta ottenute le due finestre operative, 
rispettivamente spaziale e temporale, una prima fase di 
stima del rumore 23 Snoise_est effettua una stima di un 
parametro statistic© NL n (x,y) rappresentativo del livello 
15 di rumore presente sul pixel p n (x,y) e sulla rispettiva 
finestra operativa spaziale SW n (x,y) . Di seguito ci 
riferiremo a tale rumore con il termine di rumore spaziale 
e a tale fase con il nome di stima spaziale del rumore. 

Piu in dettaglio, la prima fase di stima 23 
20 Snoise_est in primo passo effettua una stima preliminare 
di rumore (cioe di un parametro statistico rappresentativo 
del livello di rumore) calcolata localmente, cioe 
calcolata su pixel della finestra operativa spaziale 
SW n (x,y) . In un secondo passo ottiene la stima spaziale 
25 definitiva NL n (x,y) , modificando la stima preliminare in 
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base ad una stima spaziale del rumore specifica per il 
colore del pixel da filtrare e specifica per l'immagine 
Img n . 

Piu in dettaglio, la stima spaziale del rumore e 

5 ottenuta con un calcolo di tipo ricorsivo, essendo 

effettuata tenendo conto sia di un calcolo di rumore 

preliminare e locale, sia della stima spaziale del livello 

di rumore effettuata per 1' ultimo pixel filtrato 

dell' immagine Img n , avente lo stesso colore del pixel da 

10 filtrare p n (x,y) . 

In formule, nel caso in cui, ad esempio, p n (x,y) sia 

un pixel verde, si ha: 

NL n (x,y) = NL G „(x 9 y) = (i) 
= Mx,j)xiV[^(x,j)] + (l-t„(xj))xiVi:(/;/) 

in cui l'apice "G" indica che ci si riferisce al 
15 colore verde, N[SW n (x 9 y)] e la stima preliminare calcolata 
sulla finestra operativa spaziale SW n (x,y) , k n (x 9 y) e un 
fattore moltiplicativo compreso tra zero ed uno e 
determina la forza del filtro spaziale, NL G (pp G ) e la 
stima spaziale di rumore effettuata per il pixel verde pp G 
20 dell ' immagine Img n che immediatamente precede il pixel 
verde da filtrare p n (x,y) nell'ordine di scansione 
dell ' immagine Img n . 

Owiamente, se il pixel da filtrare p n (x,y) & H 
primo pixel del rispettivo colore ad essere filtrato 
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nell' immagine Img nt si ha a disposizione solo la stima 
preliminare N[SW n (x, y)] . In tal caso, si pone ad esempio, 
limitatamente al pixel p n (x,y) , k n (x,y) =1, oppure si 
assegna alia quantita NL G (pp G ) un valore arbitrario, 
5 pref eribilmente piccolo. 

La prima fase di stima spaziale del rumore 23 
Snoise_est # pud essere ad esempio eseguita come descritto 
dettagliatamente nella domanda di brevetto europeo No. 
0183 0562.3 della Richiedente, qui incorporata nella sua 

10 interezza per effetto di riferimento. Da tale documento 
risulta evidente anche il calcolo ed il significato delle 
quantita: NL n (x,y) , N[SW n (x, y)\ , k n (x,y) , NL G n (pp G ) . Si vedano 
a tal proposito le formule (1) , (2), (3), (4), (5), (6), 
(7), (8), le figure 7, 8A, BB, 10 (e la relativa 

15 descrizione) della sopra citata domanda di brevetto No. 
01830562.3. Per un ulteriore approf ondimento sulla stima 
del rumore spaziale sopra descritta, in particolare cosi 
come espressa riell' equazione (1), si faccia riferimento 
anche al brevetto statunitense US 6,108,455. 

20 La stima spaziale del rumore NL n (x,y) cosi calcolata 

interviene nella regolazione , del grado, o forza, del 
filtraggio nel caso in cui il filtraggio del pixel p n (x,y) 
sia di tipo esclusivamente spaziale. 

Terminata la fase di stima spaziale del Snoise_est 23 

25 del rumore, una fase successiva Text_est di stima di 
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disomogeneita associa un indice di disomogeneita (o grado 
di "texture") T D (x,y) al pixel da filtrare p n (x,y) in base 
ad una mi sura di disomogeneita (o analogamente di 
omogeneita) effettuata sui pixel appartenenti alia 
finestra operativa spaziale SW n (x,y) . 

L' indice di disomogeneita T D (x,y) serve a decidere se 
il pixel p n (x,y) appartiene ad una regione omogenea o 
meno, al fine di stabilire se tale pixel (dunque la 
corrispondente finestra operativa spaziale) dovra 
contribuire o meno ad una stima di rumore di tipo 
spazio/temporale che sara descritta meglio in seguito. 

Infatti, una regione omogenea, e in grado di fornire 
informazione affidabile sul rumore effettivo presente 
nell' immagine poiche le fluttuazioni (cioe le differenze) 
fra pixel appartenenti ad una regione omogenea sono 
sostanzialmente attribuibili al rumore casuale. 

All'esperto del settore sono note differenti metriche 
per il calcolo di una misura di disomogeneita su un 
insieme di pixel, per questo motivo tale argomento non 
sara ulteriormente approf ondito . Fra queste, citiamo solo 
a scopo d'esempio le seguenti: differenza massima, 
dif ferenza minima, MAD (Media delle Differenza Assolute) , 
deviazione standard, estrazione di un parametro di 
distribuzione da un istogramma dei valori digitali dei 
pixel . 
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In una realizzazione particolarmente vantaggiosa dal punto 
di vista computazionale, il parametro k n (x,y) , utilizzato 
nella fase di stima del rumore spaziale Snoise_est 23, e 
calcolato in base ad una misura di 

omogeneita/disomogeneita. In particolare e ottenuto in 
base ad un calcolo di differenze fra il pixel da filtrare 
e gli altri pixel della finestra operativa spaziale. In 
tal caso la fase di stima di disomogeneita 24 Text_est pud 
essere inglobata nella fase di stima spaziale del rumore 
23 Snoise_est, calcolarido una unica misura di 
disomogeneita utile sia alia stima del rumore spaziale, 
sia ad associare al pixel da filtrare un indice di 
disomogeneita T D (x,y) . 

Attraverso una fase di confronto 25, si verifica se 
1' indice di disomogeneita T D (x,y) del pixel p n (x,y) e 
infer iore ad un prefissato valore di soglia T h (cioe il 
pixel appartiene ad una regione giudicata omogenea) . 

Se cio non accade si passa ad una fase successiva di 
valutazione del movimento Mbt_det 27. Viceversa, prima di 
effettuare la fase Mot_det di valutazione del movimento, e 
eseguita una seconda fase 2 6 L_STnoise_est di stima locale 
del rumore (cioe relativa al pixel da filtrare) che 
servira ad ottenere una stima globale del rumore (cioe 
relativa all'intera immagine sotto trattamento) . 

Nella fase di stima 26 L STnoise est viene stimato un 
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parametro, ad esempio statistico, rappresentativo del 
rumore presente localmente nella finestra operativa 
spaziale. Tale parametro e ad esempio calcolato come una 
deviazione standard locale sulla finestra operativa 

spaziale del pixel p n (x f y) o calcolato come un'altra 
analoga misura di energia. 

Data la finestra operativa spaziale SW n (x,y) relativa 
al pixel p n (x,y) , la deviazione standard locale <r^l\ pud 
essere calcolata in accordo alia seguente formula: 



in cui N e il nuraero di pixel appartenenti alia 
finestra operativa SW n (x,y) (in questo caso N=9) ed m e la 
media dei loro valori digitali . 

Come anticipato, e come sara spiegato piu 
appro f ondit amen te in seguito, una volta terminato il 
filtraggio dell ' immagine Img nt le numerose deviazioni 
standard locali a^l\ cosi calcolate per i pixel 
dell ' immagine Img n giudicati appartenenti a regioni 
omogenee saranno utilizzate per aggiornare una stima 
globale di rumore, che chiameremo anche stima di 

rumore di tipo spazio/temporale . In particolare, tale 
stima globale sara utilizzata nel filtraggio 
spazio/temporale dell' immagine successiva Img n+ i. 
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La fase di valutazione del movimento 27 Mot_det 
confronta pixel della finestra operativa temporale 
TW n (x,y) con pixel della finestra operativa spaziale 
SW n (x,y) al fine di rilevare la presenza di movimento fra 
5 le due finest re operative, ed eventualmente valutarne 
l'entita; 

Com'e noto ad un esperto del settore, in un approccio 
di tipo u non compensato in movimento" e necessario 
prendere opportune precauzioni per evitare di introdurre, 
10 con il filtraggio spazio/temporale, artefatti dovuti al 
movimento fra immagini consecutive. In particolare, 
occorre verificare che le due finestre operative non 
contengano dati incongruenti , a causa del movimento fra 
immagini consecutive o fra porzioni di queste. 
15 Ad esempio puo accadere che una finestra operativa 

contenga pixel appartenenti ad un oggetto, mentre l'altra 
finestra contenga pixel appartenenti alio sfondo poiche 
l'oggetto si e spostato tra un'immagine e la successiva.. 

Esempi tipici di artefatti che potrebbero prodursi in 
20 questi casi sono la presenza di fastidiose scie e la 
presenza della cosiddetta u immagine fantasma" visibile 
nell' immagine filtrata e dovuta ad informazione residua di 
immagini precedenti. 

In una realizzazione preferita, la fase di 
25 rilevamento del movimento Mot det calcola come mi sura di 
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movimento una misura M(x,y) di SAD (Somma di Differenze 
Assolute) fra pixel della finestra operativa temporale e 
pixel della finestra operativa spaziale. Tale misura 
M(x,y) ha un valore tanto piu elevato quanto maggiori sono 
le differenze tra le due finestre operative e puo dunque 
essere rappresentativa del movimento fra le due finestre. 

La misura di SAD e nota al tecnico del ramo e, 
pertanto, risulta evidente dalla precedente descrizione. 

In una variante realizzativa parti colarmente 
vantaggiosa, la misura di movimento M(x,y) e una misura di 
SAD "modif icata" . Per calcolare tale misura si esegue una 
differenza in valore assoluto pixel per pixel fra le due 
finestre operative ottenendo una finestra operativa 
differenza DW n (x,y): 

DW n (x,y)=\SW n (x,y) - TW n (x,y)\. (3) 

Si calcola poi la media W avg dei pixel delle finestra 
operativa differenza DW n (x,y) , e si ottiene una misura di 
SAD modi f icata M(x,y) data da: 

M (x, y) =SAD (DW n (x, y) - W avg ) m 

Questa misura di SAD, cosi modif icata, consente 
vantaggiosamente di evitare che un cambiamento delle 
condizioni d' illuminazione fra le due immagini sia 
interpretato erroneamente come un movimento. 

In un'altra variante realizzativa particolarmente 
vantaggiosa, al fine di evitare che il valore di SAD sia 
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troppo sensibile ai valori digitali dei pixel delle due 
finestre operative (che sono contaminati da rumore) si puo 
scegliere di quantizzare leggermente i valori dei pixel, 
portando ad esempio la precisione sul valore del pixel per 
il calcolo della SAD da otto a sette bit. 

Attraverso una fase di confronto 28 , si verifica se 
la misura di movimento M(x,y) del pixel p n (x,y) e 
superiore ad un prefissato valore di soglia M h . 

Se cio accade viene effettuata una fase di filtraggio 
29 S_filter di tipo esclusivamente spaziale poiche e stato 
osservato un cambiamento eccessivo fra le due finestre 
operative. 

II filtraggio esclusivamente spaziale produce il 
pixel filtrato f_p n (x,y) a partire dai pixel della 
finestra operativa spaziale. La forza del filtraggio e 
regolata dalla stima del livello di rumore spaziale 
NL n (x,y} calcolata nella fase di stima del rumore 23 
Snoise_est. In una realizzazione preferita, il filtraggio 
digitale spaziale e effettuato secondo la tecnica 
descritta nella succitata domanda di brevetto europeo 
della Richiedente No. 01830562.3, in cui il pixel filtrato 
e ottenuto come media pesata dei pixel della finestra 
operativa spaziale (si veda in particolare la formula 
(9) ) . 

Owiamente, il filtraggio esclusivamente spaziale 
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viene inoltre effettuato per tutti i pixel della prima 
irnmagine della sequenza, poiche per questa non ci ancora 
sono dati temporali a disposizione . 

Se invece il movimento misurato M(x,y) e inferiore al 
5 prefissato valore di soglia M h viene effettuata una fase 
di filtraggio digitale ST_filter 3 0 di tipo 
spazio/temporale rappresentata in maggior dettaglio in 
f igura 6 . 

In questa, una prima fase di filtraggio 33 
10 Duncan_filt, produce un pixel filtrato prowisorio 
d_p n (x,y) in accordo ad una tecnica di filtraggio nota con 
il nome di filtraggio di Duncan che verra descritta piu 
dettagliatamente in seguito. In questa fase, il pixel 
filtrato prowisorio d j? n (x,y) e ottenuto a partire da un 
15 sotto-insieme di pixel appartenenti ad entrambe le 
finestre operative. 

Come si vede in f igura 6, se il movimento rilevato e 
giudicato suf f icientemente piccolo, cioe inferiore ad un 
ulteriore valore di soglia M± prefissato, inferiore al 
20 valore di soglia M ht sul pixel filtrato prowisorio non 
viene eseguita piu alcuna operazione e si pone: 

£jp n (x,y) = d_p n (x,y) . (4) 
In questo caso dunque il pixel p n (x,y) e 
ef f ettivamente filtrato con un metodo di filtraggio 
25 spazio/temporale di Duncan. 
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Viceversa, cioe nel caso in cui il movimento rilevato 
non e trascurabile, viene eseguita una ulteriore 
elaborazione 35 Smooth_f ilt , che produce dal pixel 
filtrato prowisorio d_p n (x,y) il pixel filtrato 
5 "definitive*" f_j> n (x,y) in accordo ad una operazione di 
w ammorbidimento" cosi come e espressa dalla seguente 
formula: 

/-Afej) B A^.p a (^)+a-A) x AM < 5 > 

in cui P n e una fattore moltiplicativo, compreso tra 

10 0 ed 1, che pud dipendere dalla misura di movimento M(x,y) 
oppure pud essere uguale per tutte le immagini della 
sequenza. Ad esempio, in una realizzazione preferita si ha 
che fi tt e circa uguale a 0,75. 

A seguito dell' operazione d' ammorbidirnento, il pixel 

15 filtrato definitive si ottiene da una porzione (in questo 
esempio 75%) del valore fornito dal filtraggio di Duncan e 
da una porzione (25%) del valore del pixel non filtrato. 
In altre parole il pixel filtrato definitivo si ottiene da 
una somma fra una frazione del pixel filtrato prowisorio 

20 ed una frazione del pixel non filtrato. 

Questo accade perche, nel caso in cui ci siano 
incongruenze non trascurabili (dovute al movimento) fra le 
finestre operative, e importante poter "trascurare" in 
misura maggiore 1 ' inf ormazione temporale e dare un po' piu 

25 di rilevanza all'immagine corrente da filtrare. 
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Verra di seguito descritta un' implementazione della 
fase di filtraggio di Duncan particolarmente vantaggiosa. 

Nel filtraggio di Duncan e dapprima selezionato un 
sotto-insieme o "range" di pixel appartenenti alle due 
5 finestre operative tramite il cosiddetto "Duncan Range 
Test" o W DRT" (si veda a tal proposito "Multiple range and 
multiple f- tests", D.B . Duncan, Biometrics, vol. 11., 
pagg. 1-42 , 1955) . 

L' applicazione del "DRT" al filtraggio digitale e 
10 noto ad esempio dalla domanda di brevetto EP 1 100 260 Al 
della Richiedente. 

Le operazioni successive di filtraggio del pixel 
p n (x,y) sono poi eseguite solo con pixel che appartengono 
al "range" selezionato tramite il DRT. 
15 Lo scopo della selezione operata con il DRT e quello 

di escludere pixel, appartenenti alle finestre operative, 
il cui valore e corrotto in misura eccessiva dal rumore. 

Ad esempio, questi pixel possono essere presenti per 
effetto di un rumore, noto con il nome di "salt and pepper 
20 noise", che e in grado di portare il valore digitale di 
qualche pixel al massimo o al minimo valore della scala 
dei possibili valori digitali. 

Inoltre lo scopo della selezione con il DRT e quello 
di escludere anche pixel che sono molto differenti dal 
25 pixel da filtrare p n (x,y) i per esempio a causa di un 



27 

differente contenuto informative Si pensi ad esempio al 
caso in cui il pixel da -filtrare p n (x,y) sia relativo ad 
uno "spigolo" e invece ci siano pixel nelle finestre 
operative appartenenti alio sfondo della scena. Se i pixel 
5 dello sfondo non fossero esclusi, il filtraggio 
comporterebbe una notevole perdita di definizione 
dell ' immagine . 

Una selezione effettuata con il DRT deve individuare 
un intervallo SI (intervallo di selezione) di valori 

10 digitali avente ampiezza S opportuna che contiene il 
maggior numero di pixel simili (in questo caso 
appartenenti alle finestre operative) al pixel da filtrare 
p n (x,y) . L' ampiezza S dell' intervallo di selezione SI e 
correlata alia deviazione standard del rumore da filtrare, 

15 che si suppone essere nota. 

Non e necessario che tale intervallo sia centrato 
attorno al pixel da filtrare: se cio accadesse, e se il 
pixel da filtrare fosse, per esempio, altamente corrotto 
da rumore il test escluderebbe pixel utili al filtraggio. 

20 In una realizzazione preferita, la deviazione 

standard di rumore utilizzata per la selezione con il DRT, 
durante il filtraggio del pixel p n (x,y) nell ' immagine lmg n , 
e la stima globale di rumore spazio/temporale cr^ L 
calcolata nel filtraggio dell ' immagine precedente Img n -i. 

25 Tale scelta e particolarmente vantaggiosa dal punto di 
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vista dell ' ottimizzazione delle risorse computazionali , in 
tal modo, infatti, si evita di dover effettuare una 
scansione completa dell ' immagine da filtrare Img n (che 
servirebbe soltanto per stimarne il rumore) prima di del 
5 filtraggio vero e proprio. 

La teoria della selezione DRT fornisce tutti gli 
strumenti per calcolare a partire dalla deviazione 
standard del rumore, 1' ampiezza S e per determinare 
l'intervallo di selezione SI ottimo. Tali strumenti 

10 richiedono tuttavia un' implement azione computazionalmente 
onerosa che mal si concilia con requisiti stringenti 
imposti nell'elaborazione di immagini in tempo reale. 

In una realizzazione preferita, un ottimo compromesso 
fra affidabilita del risultato e complessita 

15 computazionale si ottiene effettuando la selezione del 
sotto-insieme di pixel con uno fra i tre seguenti 
intervalli 511,312,813, come mostrati in figura 7: 

• l'intervallo SI1 di ampiezza S centrato attorno al 
valore digitale DV del pixel da filtrare p n (x,y) ; 

20 • l'intervallo SI2 di ampiezza S centrato attorno al 

valore digitale DV=p n (x,y) + &% L ; 

• l'intervallo SI3 di ampiezza S centrato attorno al 
valore digitale DV=p n (x,y) -a* L . 

Fra i tre intervalli SI1,SI2,SI3 si sceglie quello che 
25 contiene il numero maggiore di pixel, che in figura 7 e 



29 

l'intervallo SI1. In tal modo, si hanno buoni risultati 
anche quando il pixel da filtrare p n (x,y) e un pixel molto 
rumor o so . 

Inoltre, sempre in una realizzazione preferita, il valore 
dell'ampiezza S dell ' interval lo di selezione e calcolato 
come : 

S=3xa™ . (6) 
In una variante realizzativa, piu onerosa 

computazionalmente, ma con prestazioni ottimizzate, la 
scelta dell' intervallo di selezione e effettuata 
utilizzando opportune funzioni di peso, in accordo al 
metodo descritto nella sopra citata domanda di brevetto 
europeo EP 1 100 260 Al (in questa si vedano in 
particolare le figure lb e 4) . 

Una volta individuati, con il metodo DRT, i pixel Pj 
piu simili al pixel da filtrare appartenenti alle due 
finestre operative e compresi nell ' intervallo di selezione 
SI, si calcola il pixel filtrato prowisorio d _jp n (x,y) 
come media pesata di questi, in formule: 

PjeSI 

in cui, pref eribilmente, i coefficienti di peso aj sono 
calcolati come nella sopra citata domanda di brevetto 
europeo EP 1 100 260 Al (in questa si veda in particolare 
pag. 6, righe 41-50 ). 
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Ritornando al diagramma di figura 3, dopo il 
filtraggio, che come e stato spiegato awiene in accordo 
alia fase 29 S__filter o in accordo alia fase 3 0 ST_filter, 
una fase 31 di controllo verifica se il pixel p n (x,y) 
5 appena filtrato e 1' ultimo dell ' immagine Img n . Se p n (x,y) 
non e 1' ultimo pixel, il metodo rappresentato come 
successione di fasi in figura 3 viene applicato al pixel 
successivo nell'ordine di scansione, ad esempio al pixel 
p n (x,y+l) . 

10 Altrimenti, viene eseguita una fase di stima del 

rumore globale di tipo spazio temporale sulla base delle 
numerose stime local i cr^f calcolate, per i pixel 
dell' immagine Img n giudicati appartenenti a region! 
omogenee, nella fase di stima del rumore locale 26. In 

15 particolare, tali stime sono utilizzate per aggiornare una 
stima globale <x^f, di rumore di tipo spazio/temporale che 
sara poi utilizzata nel filtraggio spazio/temporale 
dell' immagine successiva Img n+1 . Pref eribilmente la stima 
globale crf+, e rappresentativa della deviazione standard 

20 del rumore ed calcolata come media delle numerose stime 
locali (deviazioni standard) cr^ . 

In alcune situazioni pud accadere pero che immagini 
adiacenti abbiano valori di deviazione standard globali 
<t^j eccessivamente discordanti fra loro e di conseguenza 

25 il filtraggio spazio/temporale potrebbe filtrare immagini 
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adiacenti con intensita troppo differente. Questo 
genererebbe un fastidioso sfarfallio (flickering) visibile 
nella riproduzione della sequenza. 

In una realizzazione preferita, per evitare questo 
5 inconveniente, la deviazione standard globale crf^, , 
calcolata come media di deviazioni standard locali, viene 
modificata effettuando una media temporale in modo 
ricorsivo fra un certo numero (ad esempio, due) di 
deviazioni standard globali a GL relative ad immagini 
10 consecutive. In f ormule : 

*%=r*°Z+a-r)*<r? (8) 

in cui y e un numero compreso tra 0 ed 1, a^ L la stima di 
rumore globale aggiornata durante il filtraggio 
dell' immagine precedente Imgn-i ed utilizzata nel 

15 filtraggio spazio/temporale dell' immagine corrente Img n . 
Ad esempio, il numero y pud essere circa uguale a 0,75. 

Verranno ora descritte forme di realizzazione 
alternative al particolare metodo sopra descritto con 
riferimento alia figura 3 

20 In una variante realizzativa, con riferimento alia 

figura 3, il blocco di valutazione del movimento 27 
Mot__det si limita ad effettuare una valutazione di 
presenza/assenza di movimento fornendo in uscita una 
misura M(x,y) di tipo binario ed indicativa della 
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presenza/assenza di movimento. Ad esempio, si ha in uscita 
M(x,y)=l se viene rilevata la presenza del movimento, 
M(x,y)=0 altrimenti. In questo caso il valore di soglia M h 
pud essere ad esempio posto uguale a 0,5. 

La fase 24 Mot_det puo rilevare la presenza di 
movimento con un metodo di "rilevazione di scia" che opera 
sempli cement e calcolando le differenze pixel per pixel fra 
le due finestre operative. Se tali differenze hanno tutte 
lo stesso segno, posit ivo o negative, viene rilevata la 
presenza di movimento, altrimenti viene rilevata l'assenza 
di movimento . 

Nel caso in cui sia rilevata la presenza del 
movimento viene effettuata la fase 29 di filtraggio 
spaziale, gia descritta in precedenza. 

Nel caso in cui invece sia rilevata l'assenza di 
movimento, con riferimento alia figura 6, viene effettuata 
una fase di filtraggio spazio/temporale 3 0 ST_filter, in 
cui il pixel filtrato e ottenuto con una fase di 
filtraggio di Duncan 33 Duncan_Filt indif f erentemente 
seguita, o meno, da una fase di filtraggio di 
ammorbidimento 3 5 Smooth_Filt. 

In un'ulteriore variante realizzativa le due fasi di 
filtraggio S_filter 29 e ST__filter 30, rispettivamente 
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spaziale e spazio temporale, effettuano le medie pesate 
per il calcolo del pixel filtrato f_p n (x,y) sostituendo ad 
alcuni pixel che intervengono nella media pesata il 
rispettivo valore filtrato, nel caso in cui questo si gia 
5 disponibile nel buffer d'uscita. 

Risultati sperimentali hanno mostrato che il metodo 
di filtraggio proposto e in grado di fornire concreti 
vantaggi in termini di qualita d'immagine e in termini di 
efficienza di codif ica/conrpressione, al contempo 

10 richiedendo un' allocazione modesta di risorse 
computazionali e di memorizzazione. Questo rende 
l'utilizzo del metodo della presente invenzione 
particolarmente vantaggioso in applicazioni che richiedono 
capacita d' elaborazione in tempo reale. 

15 In figura 8, la curva indicata con n CFA_f iltered" 

riporta, immagine per immagine, una misura di qualita per 
una sequenza di trecento immagini filtrate secondo la 
presente invenzione. La misura e effettuata su immagini 
filtrate ed interpolate. 

20 Sempre in figura 8, la curva indicata con "Noisy" 

riporta la stessa misura effettuata sulla stessa sequenza 
non filtrata; anche in questo caso la misura e effettuata 
su immagini interpolate . 
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I risultati di figura 8 si riferiscono in particolare 
ad una misura di PSNR (Peak-toPeak Signal to Noise Ratio) . 
II PSNR e una misura di tipo standard ed e rappresentativa 
della qualita di una immagine, in particolare indica la 
quantita di segnale presente in una immagine rispetto alia 
quantita di rumore . 

Dai risultati sperimentali di figura 8 si osserva che 
la sequenza filtrata presenta una misura di PSNR piu 
elevata (tipicamente il guadagno £ di circa 3 dB) , il che 
e indicativo di una qualita maggiore. 

In figura 9 sono mostrati due schemi di elaborazione 
per produrre a part ire da una sequenza di immagini CFA 
rumorose una sequenza di immagini filtrate e compresse / 
codificate in accordo alio standard MPEG4 ♦ 

In particolare, nel cammino superiore, indicato con 
El, le immagini CFA sono filtrate con un metodo di 
filtraggio secondo la presente invenzione dal filtro 5 CFA 
NF, successivamente sono interpolate dal blocco 7 IGP e 
codificate / compresse dal blocco 8 MPEG4 - Encoder . 
Nel cammino inferiore, indicato con E2, le immagini CFA 
rumorose sono prima interpolate dal blocco 7 IGP, poi 
secondo uno schema di filtraggio convenzionale, sono 
filtrate con un metodo di filtraggio compensato in 
movimento effettuato immediatamente prima della codifica 
compressione MPEG ed indicato con MC_NF . In particolare il 
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metodo di filtraggio MC_NF utilizza la stima di movimento 
e la compensazione del movimento del codificatore MPEG per 
effettuare un filtraggio digitale di tipo spazio / 
temporale compensato in movimento. 
5 Risultati sperimentali hanno dimostrato che lo schema 

d'elaborazione El, dunque in accordo alia presente 
invenzione, consente di ottenere un guadagno del 20% in 
termini di u bit rate" rispetto alio schema di elaborazione 
conventionale E2 . 

10 Questo pud essere spiegato considerando il fatto che 

nello schema d'elaborazione E2, la stima del movimento e 
la compensazione del movimento awengono su immagini 
affette da rumore, dunque possono non fornire un risultato 
ottimale, pregiudicando cosi sia la qualita delle immagini 

15 che 1'efficienza d compressione . 

II metodo di filtraggio secondo 1 ' invenzione, sopra 
descritto in una sua forma di realizzazione preferita, pud 
essere implementato utilizzando hardware, software, o una 
combinazione di hardware e software. In quest 'ultimo caso 

20 il metodo pud essere implementato in uno specifico 
circuito integrate ASIC (Application Specific Integrated 
Circuit) . 

Nel caso in cui sia implementato in un dispositivo 
per acquisire sequenze di immagini, il metodo secondo la 
25 presente invenzione pud vantaggiosamente essere eseguito 
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mediante risorse d' elaborazione (ad esempio, DSP) 
condivise con altre applicazioni all'interno dello stesso 
dispositivo. 

Owiamente, al metodo sopra descritto un tecnico del 
ramo, alio scopo di soddisfare esigenze contingenti e 
specif iche, potra apportare numerose modifiche e varianti, 
tutte peraltro contenute nell'ambito di protezione 
dell ' invenzione, quale definito dalle seguenti 
rivendicazioni . 
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RIVENDICAZIONI 

1. Metodo per filtrare il rumore di una sequenza di 
immagini digitali ( Img 1/W , Img n , Img B +i,—) in formato CFA 
comprendente le seguenti fasi: 

- processare una prima immagine CFA {Img n - ± ) della 
sequenza per ottenere una corrispondente immagine CFA 
migliorata {£_Img n -i) avente rumore ridotto; 

- elaborare almeno un pixel (p n (x,y)) di una seconda 
immagine CFA (Img n ) della sequenza che segue temporalmente 
detta prima immagine (X/ngn-j) , detta fase di elaborare 
1 ' almeno un pixel fornendo un corrispondente pixel 
filtrato (f_p n (x,y)) ed includendo le seguenti operazioni : 

selezionare (SW_sel) un primo insieme di pixel 
(SW n (x, y) ) comprendente detto almeno un pixel ed una 
pluralita di pixel della seconda immagine CFA ad esso 
spazialmente adiacenti; 

- selezionare (TW_sel) un secondo insieme di pixel 
(TW n (x,y)) comprendente pixel dell ' immagine CFA migliorata 
omologhi ai pixel di detto primo insieme; 

- effettuare un filtraggio digitale di un primo tipo 
(ST_filter) utilizzando pixel appartenenti a detti primo e 
secondo insieme per generare il pixel filtrato 
(f_j>n (x,y) ) ■ 

2. Metodo secondo la rivendicazione 1, comprendente 
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inoltre una f ase di : 

effettuare (Mot_det) una prima valutazione di 
movimento (M(x,y)) dell'almeno un pixel, utilizzando pixel 
appartenenti a detti primo e secondo insieme; 

ed in cui detto almeno un pixel e tale che detta 
valutazione di movimento e inferiore ad un primo valore di 
soglia (M h ) . 

3. Metodo secondo la rivendicazione 1 comprendente 
inoltre una fase di: 

- elaborare un ulteriore pixel {p n (x,y)) di detta 
seconda immagine, detta fase di elaborare 1' ulteriore 
pixel fornendo un corrispondente ulteriore pixel filtrato 
ed includendo le seguenti operazioni: 

selezionare (SW_sel) un terzo insieme di pixel 
comprendente detto ulteriore pixel ed una pluralita di 
pixel della seconda immagine CFA ad esso spazialmente 
adiacenti; 

- selezionare (TW_sel) un quarto insieme di pixel 
comprendente pixel dell' immagine CFA migliorata omologhi 
ai pixel di detto terzo insieme; 

- effettuare (Mot_det) una ulteriore valutazione di 
movimento {M(x,y)) dell ' ulteriore pixel, utilizzando pixel 
appartenenti a detti terzo e quarto insieme; 

se detta ulteriore valutazione di movimento e 
inferiore a detto primo valore di soglia (M h ) , effettuare 
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un filtraggio digitale di un secondo tipo (S_filter) che 
genera l'ulteriore pixel filtrato utilizzando pixel 
esclusivameiite appartenenti a detto terzo insieme. 

4. Metodo secondo la rivendicazione 1, in cui ciascuna 
5 inunagine CFA della sequenza e formata da una rispettiva 

matrice di pixel, i pixel di detta matrice essendo 
associati, in base alia rispettiva posizione nella 
matrice, ad una di un insieme di componenti cromatiche ed 
in cui detto primo e detto secondo insieme comprendono 
10 pixel associati alia stessa componente croraatiea 
dell'almeno un pixel. 

5. Metodo secondo la rivendicazione 4, in cui le 
immagini della sequenza sono in formato CFA di Bayer, e 
dette componenti cromatiche appartengono all 'insieme 

15 comprendente il colore rosso, il colore verde ed il colore 
blu. 

6. Metodo secondo la rivendicazione 5, in cui la fase di 
selezionare il primo insieme di pixel e effettuata 
mediante matrici di selezione (SM_G, SM_R, SM_B) differenti 

20 a seconda della componente cromatica dell'almeno un pixel, 
le matrici di selezione essendo tali da selezionare pixel 
disposti in un intorno dell'almeno un pixel ed aventi la 
stessa componente cromatica di questo e scartando gli 
altri, le matrici di selezione essendo inoltre identiche 

25 per le componenti cromatiche rosse e blu. 
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7. Metodo secondo la rivendicazione 1, inoltre 
includente una fase (G_STnoise_calc) di effettuare una 
stima di un parametro statistico cr™ rappresentativo del 
rumore globale presente in detta prima immagine ( Imgn-x) , 

5 il filtraggio digitale del primo tipo utilizzando inoltre 
detto parametro statistico. 

8. Metodo secondo la rivendicazione 7, comprendente 
inoltre le seguenti . f asi : 

selezionare una pluralita di pixel della prima 
10 immagine; 

- calcolare una pluralita di stime locali, calcolando 
per ciascun dato pixel di detta pluralita di pixel una 
rispettiva stima di un parametro statistico 
rappresentativo del rumore locale presente in un intorno 
15 del dato pixel; 

ed in cui detta stima del parametro statistico del rumore 
globale cr£ L e ottenuta a partire da detta pluralita di 
stime locali . 

9. Metodo secondo la rivendicazione 8, in cui dette 
20 stime locali sono misure locali di varianza. 

10. Metodo secondo la rivendicazione 8, in cui detta 
pluralita di pixel comprende pixel appartenenti a regioni 
omogenee della prima immagine. 

11. Metodo secondo la rivendicazione 1, il cui detto 
25 filtraggio digitale del primo tipo utilizza un sotto- 
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insieme di pixel appartenenti a detto primo e da detto 
secondo insieme di pixel, detto sotto- insieme essendo 
individuato mediante una ulteriore fase di selezione 
effettuata in accordo al Duncan Range Test. 
5 12. Metodo secondo la rivendicazione 3, comprendente 
inoltre una fase di stimare per 1' ulteriore pixel un 
ulteriore parametro statistico rappresentativo del rumore 
locale presente sui pixel di. detto . terzo insieme,. detto 
ulteriore parametro essendo stimato anche in base alio 
10 specif ico colore dell ' ulteriore pixel, il filtraggio 
digitale del secondo tipo utilizzando detto ulteriore 
parametro. 

13. Filtro, per ridurre il rumore di una sequenza di 
immagini in formato CFA, caratterizzato dal fatto detta 

15 riduzione del rumore e effettuata secondo il metodo in 
accordo ad una qualsiasi delle precedenti rivendicazioni . 

14. Dispositivo di acquisizione per acquisire una 
sequenza di immagini digitali in formato CFA, comprendente 
un sensore per acquisire dette immagini, detto sensore 

20 comprendendo un filtro CFA, caratterizzato dal fatto che 
la sequenza di immagini in formato CFA e elaborata secondo 
il metodo di filtraggio di rumore in accordo ad una 
qualsiasi delle rivendicazioni da 1 a 12 . 
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M Metodo per filtrare il rumore di una sequenza di 
immagini digitali" 
RIASSUNTO 

Metodo per filtrare il rumore di una sequenza di immagini 
digitali Img lf Img n , Img n+1 ,„. in formato CFA comprendente le 
seguent i f as i : 

processare una prima immagine CFA img^x della 
sequenza per ottenere una corrispondente immagine CFA 
migliorata f^Imgn-x avente rumore ridotto; 

- elaborare almeno un pixel p n (x,y) di una seconda 
immagine CFA lmg n della sequenza che segue temporalmente 
detta prima immagine lmgn-i, detta fase di elaborare 
1' almeno un pixel fornendo un corrispondente pixel 
filtrato f_p n (x,y) ed includendo le seguenti operazioni: 

selezionare SW_sel un primo insieme di pixel 
SW n (x,y) comprendente detto almeno un pixel ed una 
pluralita di pixel della seconda immagine CFA ad esso 
spazialmente adiacenti ; 

selezionare TW_sel un secondo insieme di pixel 
TW n (x,y) comprendente pixel dell' immagine CFA migliorata 
omologhi ai pixel di detto primo insieme; 

effettuare un filtraggio digitale di un primo tipo 
ST__filter utilizzando pixel appartenenti a detti primo e 
secondo insieme per generare il pixel filtrato f _Pn(x r y) 
(Figura 3) 
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